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Abstract: 

We have computed the hght quark masses using the 0(a^) improved Alpha action, in the 
quenched approximation. The renormahzed masses have been obtained non-perturbatively. 
By eliminating the systematic error coming from the truncation of the perturbative series, 
our procedure removes the discrepancies, observed in previous calculations, between the 
results obtained using the vector and the axial-vector Ward identities. It also gives values 
of the quark masses larger than those obtained by computing the renormalization constants 
using (boosted) perturbation theory. Our main results, in the RI (MOM) scheme and at 
a renormalization scale /i = 2 GeV, are mf^ = 138(15) MeV and mf^ = 5.6(5) MeV, 
where mf^ is the mass of the strange quark and mf^ = (m^^ -|- m^^)/2 the average mass 
of the up-down quarks. From these results, which have been obtained non-perturbatively, 
by using continuum perturbation theory we derive the MS masses, at the same scale, and 
the renormalization group invariant (m^*^^) masses. We find tt^^-^o m5 _ 121(13) MeV 
and mf^^^ = 4.9(4) MeV at the next-to-leading order; mf^^^ = 111(12) MeV, 
^N^LOMS ^ 4 5(^4) j^gY^ ^RGi ^ 177(19) MeV and mf^^ = 7.2(6) MeV at the next-to- 
next-to-leading order. 
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1 Introduction 



The values of the quark masses obtained from lattice simulations have recently attracted 
the attention in the physics community. The reason is that, for these quantities, the 
lattice method is unique: from the non-perturbative calculation of hadronic quantities, 
this approach allows in fact a consistent determination of the quark masses, defined as 
effective couplings renormalized at short distances. Yet, with the errors quoted by the 
authors the values of the light and strange quark masses computed in different 

numerical simulations (mainly in the quenched approximation) are often in disagreement. 
The differences originate mainly from two sources: on the one hand, from the different 
procedures used to compute the renormalized mass from the bare lattice one; on the other, 
from the different methods used to extrapolate results, obtained at finite values of a, to 
the continuum limit. 

Our numerical calculations have been done with the non-perturbatively improved fer- 
mion action (which we will denote as the Alpha action) introduced in ref. [^, see also 
T^-|T^. The use of a non-perturbatively improved action and operators reduces dis- 



cretization errors to (9(a^). Thus, at least for light quarks (namely the u, d and s quarks), 
we expect these errors to be rather small, i.e. much smaller than other systematic ef- 
fects. The other source of uncertainty arises from the truncation of the perturbative series 
in the definition of the renormalized mass. This problem can be eliminated by using a 
non-perturbative method for renormalizing the lattice operators. In our study, we have 
computed the relevant renormalization constants Zg and Zp, in the chiral limit, using the 
non-perturbative renormalization procedure on quark states proposed in ref. (in the 
following we denote this method as NPM). As a check of the accuracy of the NPM, we have 
also computed the renormalization constants of the vector and axial-vector currents and 



compared the results with those obtained in refs. |]T4|, see also [T^. At /? = 6.0 and 6.2, 
our results for these currents agree within less than 5% with previous determinations. This 
makes us confident that the values of Zs and Zp determined with the NPM are correct. 
Note that, at the same values of /3, the value of Zp {Zs) computed in boosted perturbation 
theory is larger than our estimate by about 30% (by about 10%) Q This implies that the 
use of perturbation theory leads to an underestimate of the quark mass by an amount much 
larger than the expected discretization errors 0. It also makes the two determinations of 
the masses, from the vector and axial vector Ward identities, different. With the NPM, 
instead, we find essentially the same value of the quark masses with the two methods. 

This letter is focused on the calculation of the light quark masses: we describe the 
procedure followed in the determination of the bare masses, m, discuss the definition of 
the renormalized and of the renormalization group invariant (RGI) ones, rh and m^*^^, 
and present the final results and errors for these quantities. All details of the (standard) 
analysis of the light hadron spectrum and decay constants, together with a study of the 



^ We used for the boosted coupling ||T^, g'^ = gQ/{P), where the average plaquette at /? = 6.2 
(/3 = 6.0) is given (P) — 0.6136 {{P) — 0.5937) as inferred by our simulations. 

^ The observation that the NPM gives larger quark masses was first made in ref. , and then confirmed 
in a systematic study at different values of /? and with different actions in ref. j7| . 
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energy-momentum relation, can be found in ref. For the light hadron spectrum, the 

calibration of the lattice spacing, the values of the Wilson parameter in the chiral limit, 
Hcrit, and the values of k corresponding to the light (up and down) and strange quark 
masses, we substantially agree with previous studies of the same quantities We 
performed two independent simulations: a test run at (3 = 6.0 on a 16'^ x 32 lattice (with a 
sample of 45 configurations at four values of the light quark masses) and a run at = 6.2 
on a 24^ x 64 lattice (on a sample of 100 configurations at four values of the light quark 
masses). The final results are based on the "data" at /3 = 6.2. Those of the test run have 
only been used as a consistency check of the stability of the physical values of the quark 
masses. We will compare the results at (3 = 6.0 with those at /? = 6.2 at the end of the 
paper. 



2 Definition of the quark masses 

We now explain the procedure followed to extract the physical values of the quark masses. 
The starting point are the vector and axial-vector Ward identities (WI) ^ 

V,{aV,\P) = im,{fx)-m2ifx)){a\Sifxm , (2.1) 

V^(a|4l/3) = (mi(/i) +m2(/i))(a|P(^)|/5) , (2.2) 

where |a) and represent generic external physical states; and are the normalized 
currents which obey the current algebra commutation relations; S{fi), P{fi) and rhi{fi) 
are operators and quark masses renormalized at the scale /i in a given scheme. Note 
that the products rh{fi)S{fi) and m{fi)P{fi) are regularization, renormalization and scale 
independent. 

Ward identities as (|2.1| ) and ( p.2|) can also be written for the lattice currents. In this case 
they are valid up to terms of C(a), 0{asa) or O(a^), depending whether we use the Wilson, 
the tree-level improved or non-perturbatively improved actions and operators 0. 

In order to determine the renormalized masses it is obvious that we have first to specify 
the scheme used to renormalize the scalar and pseudoscalar densities. Given the uncer- 
tainties of lattice perturbation theory, we have renormalized these quantities by imposing. 



non-perturbatively, the following renormalization conditions 

{q{p)\S{^^Mp)) ^ Zs{fi){q{p)\S\q{p)) = 1, 

{q{p)\P{^iMp)) ^ Zpif^){q{p)\P\Q{p)) = 1, (2.3) 

in the chiral limit. In the above equations \q{p)) represents an external off-shell quark 
state with virtuality = /i^, and the matrix elements are evaluated in the Landau gauge. 
The conditions (|2.3|) at large values of fi"^ ensure that S and P belong to the same chiral 
multiplet, i.e. they obey to the relevant Ward identities. It can be shown that this is true 

On the lattice is written in terms of (improved) finite differences. 
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also at the improved level, i.e. that, using the improved action of ref. the matrix 
elements of S and P renormalized as in eq. ( p.3| ) have discretization errors of O(a^) • 

Since vector symmetries are preserved by the lattice regularization, unlike axial vector 
ones, the procedures which are usually employed to extract the quark masses using the 
vector or axial vector Ward identities are different and we discuss them separately. 

• Determination of quark masses from the vector Ward identity. 

/^From the lattice version of eq. one finds 



m(/i) = Zg^{fi,ma,0) 



1 /I 



1 



2(1 \ K t^crit 



(2.4) 



where, in the definition of Zg-, we added those terms which are necessary to improve 
the scalar density out of the chiral limit 



Zs{^i, mia, m2a) = Zs{fi)(l + — (mia + 



By defining Zs{fi) 
relation 



(2.5) 



Z^i^jji) and hs = —2bm, from eq. (p.4|) we derive the standard 



m(/i) = Z^(/i) m = ZmifJ') 



1 /I 1 



1 + hra— [ - - 



1 



2 a \k l^crit 



(2.6) 



We have used the above equation to determine the quark mass, with Z„ 
Zs{fi)~^ as computed from eq. (|2.3|) , Hcrit as fixed from the squared pseudoscalar 
meson mass, the calibration of the lattice spacing either from rrip or from rriK* and 
1/k = l/nu,d,s from the meson spectroscopy. The discussion of the values of Zs and 
Zp and the relative errors can be found in sec. |3|. 



As far concerned, the method of ref. [E3] has difficulties in computing it. 

The gauge invariant procedure of ref. has not been applied to date. A different 
approach to fix bm non-perturbatively has been proposed in ref. We have tried 
the same technique and found that the results for bm (and for all the analogous 
quantities such as b^, by etc.) are very unstable and that these constants cannot 
be determined reliably [^] 0. For this reason, we preferred to take bm from boosted 
perturbation theory, using 6^ = b^^+b^^g"^ = -1/2-0.0962 = -0.652 [g? 
uncertainty on the perturbative value of bm were as large as b^g"^ 
a relative error of 0{6 ■ 10~'^) and C(4 • 10"'^) on {rhd + fhu)/2 and rfis, respectively. 
Given the other sources of errors (due to the uncertainty in the determination of Zg 
and to the quenched approximation), this error is negligible and will be ignored in 
the following. 



If the 
this would induce 



In ref. 0, they call HS (from Hadron Spectrum) the method based on eq. ( p.6| ). This 
is rather misleading because the possibility of relating l/2a{l/n — l/hZcrit), or any 



''A new method for a non-perturbative determination of the 6s, based on the lattice chiral Ward iden- 
tities, has been recently presented at the Lattice '98 conference |2^. However, this method has not yet 
been implemented with the non-perturbatively improved Alpha action. 
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improved version of it, to the quark mass, entirely relies on the validity of the lattice 
vector WI. The hadron spectrum enters in this case, as in the case of the axial-vector 
WI, only because it is needed to fix the bare quark mass from the physical value of 
some hadronic quantity. 

• Determination of quark masses from the axial-vector Ward identity. 

In this case, the simplest procedure is to use the axial Ward identity computed on 
hadron states at rest, for degenerate quark masses. One gets 

= —^p - («| Vq^o + acAVo-P|/^) 

Zp{fi) Zp{fj,,ma,ma) 2a{a\P\l3) 

Zp{fi) 

where and bp play the same role for the axial current and the pseudoscalar density 
as 65 in the case of the scalar operator. In the absence of a precise non-perturbative 
determination of these quantities, we have used, as we did above for b^, boosted 
perturbation theory, namely bA = 1 + 0.15225(2 = 1.240 and bp = l + O.lSSlgf^ = 1.241 
[p^ . The same considerations made for bm about the systematic error induced by 
perturbation theory in the calculation of bA and bp on the masses of the light quarks, 
remain true in this case. We extracted the quark mass by computing the ratio in 
eq. (|2.7| ) at values of k corresponding to the physical meson masses (m^, rriK and 
m^), taking the non-perturbative determinations of Za and Zp{fi), evaluated in the 
chiral limit. 



1 + bAtna (a| Vo^o + clcaV^PI 



1 + bpTua 



2a{a\P\ 



(2.7) 



3 Non-perturbatively renormalized quark masses 

In this section, we first discuss the values and errors of the renormalization constants Zs{p), 
Zp{fj,) and Za- We then explain the procedure followed to extract the bare masses, fh and 
p, using eqs. (|2.6|) and (|2.71 ). /^From the Zs and from the values of the bare masses, we 
have obtained the values of the non-perturbatively renormalized quark masses in the RI 
scheme. Using perturbation theory, we then get the MS and RGI masses at the NLO and 
N^LO. 

The relevant Green functions (vertices and propagators) necessary to the determination 
of the Zs with the NPM have been computed with external quark states of virtuality p on 
the same configurations, and for the same values of k, as all the other quantities discussed 



in this paper. We used quark propagators improved following the strategy of refs. |2^, p6 
The Green functions have then been extrapolated (linearly in the quark masses) to the 
chiral limit. 



As explained in ref. |]T6|, the NPM is expected to work when p satisfies the condition 
^QCD ^ p ^ 1/ a. In this region, Za and Zp/Zs should be independent of /z, since the 
non-perturbative method is equivalent to the Ward identities. We monitored the behaviour 
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Figure 1: Zs obtained by using the NPM as a function of the scale fi^a"^. The dashed and 
solid curves represent the solutions of the renormalization group equations at the leading 
and next-to-leading order, respectively. 



of Za and Zp/Zs as a function of yU to find the range of yU where these quantities exhibit 
a plateau. Ai (3 = 6.2, we find that Z^ is essentially a constant for /x^a^ < 2, while Zp/Zs 
has a plateau for 1 < yU^a^ < 2, corresponding to 2.8 GeV < /i < 4.0 GeV. /,From the 
analysis of the plateaus of Za and Zp/Zs, we have obtained Za = 0.793(5), in agreement 
with the value Za = 0.809 of ref. |]14[, and Zp/Zs = 0.78(1). The latter in disagreement 
with the boosted result Zp/Zs = 0.95. Similar results, although with larger statistical 
fiuctuations, and systematic uncertainties, have been found at /5 = 6.0 p6| . 

As a consistency check of our results, we have also studied the /i dependence of Zs- In 
fig. ID we show the behaviour of Zs as a function of fi compared with the solution of the 
renormalization group equations 



Zsifi) = Zsifio) 



' ajjiy 
as(/io), 



7o/2/3o 



An 



J 



(3. 



where, in Landau Rl and with Uf = 



7o 



71 = -252 , Po 



11 



J 



102 



(3.9) 



By fitting the numerical results to eq. 



^), with Aqcd as free parameter, we get Aqcd = 

° 251 ±21 



200 ± 55 MeV in good agreement with the determination of ref. ||31 
MeV. It is reassuring that eq. 



is]) fits rather well the behaviour of Zs{jj) down to rather 
low values of iJ?. From the fit, and using Zp/Zs = 0.78(1), we find the results given in 
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(3 = 6.2 




Zs (NPM) 


Zs (BPT) 


Zp (NPM) 


Zp (BPT) 


2.80 GeV 


0.60(1) 


0.66 


0.47(1) 


0.62 


2.00 GeV 


0.55(1) 


0.61 


0.43(1) 


0.57 


(3 = 6.0 


/i 


Zs (NPM) 


Zs (BPT) 


Zp (NPM) 


Zp (BPT) 


2.00 GeV 


0.55(3) 


0.63 


0.39(3) 


0.59 



Table 1: Zs and Zp obtained by using the NPM and boosted perturbation theory (BPT), 
at two typical reference scales, /i = 2.80 GeV, corresponding to fi'^a'^ = 1 at (3 = 6.2, 
and fi = 2 GeV which is the reference renormalization scale for the quark masses in 
lattice calculations. We also used the values of Za as obtained with the NPM, namely 
Za = 0.793(5) and Za = 0.78(1) at [3 = 6.2 and 6.0 respectively. 



table 1. We stress that the use of perturbation theory induces an error of 10-30% in the 
extraction of the values of the quark masses, depending whether one uses the vector or 
axial- vector Ward identity (i.e. Zs or Zp). For completeness we also give the results at 
(3 = 6.0. More details on the analysis of the Zs obtained with the NPM can be found in 
ref. m ^ 

Eqs. ( |2.6| ) and ( p.7| ) determine the bare masses rh and p as a function of the hopping 
parameter k. The physical values of etc., are fixed, together with the lattice 

spacing, using a certain number of physical conditions. We now explain the procedure 
followed in our analysis. 

We fit the mass of the vector meson to the expression 

aMv = C + L{aMpsf + Q{aMpsY , (3.10) 

where Mps is the pseudoscalar mass. The cases with Q 7^ or Q = are denoted as 
quadratic or linear fit respectively. By taking Mps = m.„ {Mps = rriK) and My = nip 
{My = rriK*), where and nip are the experimental numbers, we determine the value 
of the lattice spacing. This determination, being based on the physical spectrum without 
reference to any definition of the quark masses, is valid up to 0{a'^). It has been called 
"lattice-plane method" (p-7r or K*-K) in ref. [0. To determine m and p, we must study 
another physical quantity. The best (and most popular) choice is Mj.^, since it vanishes 
in the chiral limit and for this reason is very sensitive to the precise value of the quark 
masses. We fit Mp^ to 

a^Mpg = Lpsa{nii + 1712) + Qpsa'^i'mi + "^2)^ , (3-11) 

'''The amputated correlation function of the pseudoscalar density, between external off-shell quark states, 
is expected to receive in the chiral limit a non-perturbative contribution from the Goldstone pole . We 
have eliminated such a contribution and determined the renormalization constant Zp by evaluating the 
ratio Zp/Zs at large values of fi^ and then calculating Zp as (Zp/Zs) ■ Zs- 
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(3.12) 



depending whether we use the vector or axial-vector WI. In the above equations, with the 
values of a fixed from eq. ( p.lO|) and by imposing Mps = 171^^ or rriK, we determine rhi {pi) 0. 



For consistency, quadratic or linear fits should be used both in eq. (|3.1CI|) and eq. ( |3.11D or 



( p.l2| ). We have also determined the strange quark mass from the vector meson mass, by 
fitting My to 

aMy = Cy + Lya{mi + 1712) + Qyo^irhi + 1712^ , (3.13) 
aMy = Cy + L'ya{pi + P2) + QyO^i^pl + p2f , (3.14) 

and using the experimental value of m^. 

^From the analysis of the results obtained using quadratic and linear fits, we reached 
the following conclusions 

• As already observed in ref. [§], the coefficient Qps is large and positive, Qps = 
1.54(14). It is very difficult to explain the sign of Qps as a consequence of the use of 
perturbation theory for such an explanation would require the perturbative value 
of bm to be wrong both in size and in sign. This is very unlikely (and to our knowledge 
never happened to be the case). Thus, we believe that the positive curvature Qps is 
a physical effect (at least in the quenched approximation). 

• In the range of masses used in our calculation, corresponding to the values of k = 
0.1333, 0.1344, 0.1349 and 0.1352 (from heavier to lighter), the quadratic fits on 
all the four quarks give consistent results and statistical errors for the quark masses. 
Moreover, the value of HcriUc is in agreement with that obtained form the axial-vector 
WI. 

• With our data, instead, a linear fit of the pseudoscalar mass on all the four values of 
K gives a value of HcriUc incompatible with that obtained from the axial-vector WI. 
This confirms that we need a quadratic fit for the pseudoscalar mass. 

• Within the accuracy of our results, the coefficient Q = —2.2 ±1.5 {Qy = —2.1 ± 2.5) 
is compatible with zero. Nevertheless, by using the linear or quadratic fit to My, the 
central value of the lattice spacing, and the value of the strange quark mass, vary by 
about 10%. 

• If we assume that the quadratic correction to My is negligible, i.e. if we take Q = 
but fit quadratically Mpg (as done in ref. |^), we obtain = 2.59(12) GeV (using 
nip-m.,,) and arhs = 0.031(3) in excellent agreement with refs. p, |^. 

• With a quadratic fit for both My and Mpg, however, we obtain a value of the 
inverse lattice spacing which is sensibly higher than 2.59 GeV, i.e. = 2.84(32). 



^ Since we always work with degenerate quarks, we are unable to fit a term oc (mi — 7712)^. In order to 
separate m; from (using and mn), we must, however, distinguish the two quark masses. 
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Although the two results look compatible within the errors, we believe that the shift 
is a systematic effect, since the two values of are determined on the same set of 
configurations. Correspondingly, the values of the quark masses are also modified. 
In the quadratic case, we get arhg = 0.026(5). We stress again that the shift of arhg 
toward smaller values, although compatible with 0.031(3), is a systematic effect. 

It is well known, see for example refs. ^, 0, that the strange quark mass extracted 
using rriK differs from that obtained from by about 15%. Moreover, the values of 
the inverse lattice spacing determined from mp-niTr or mK^mK* are slightly different 
(by about 100 MeV). These effects were found irrespectively of the action (Wilson, 
tree-level improved or non-perturbatively improved) used in the numerical calcula- 
tions and their origin can be traced from the fact that the slope L in eq. ( |3.10|) is 



smaller than its experimental value. With our data, we find the same effects if we use 
the linear fit for My (with a difference of about 10-15% for am^, for example). With 
quadratic fits, we find two nice features: on the one hand the difference in the value 
of from mp-irtT^ and mK*^mK is strongly reduced; on the other, we get about the 
same strange quark mass from rriK and m^. 

Unfortunately, within our precision, we are unable to fix Q well enough. Nevertheless, 
at the price of increasing the statistical error, we will take as best estimates of the 
quark masses those obtained by using quadratic fits for both My and Mp^. Since 
this is at present one of the largest sources of uncertainty, we believe that a precise 
measurement of the quadratic term in the dependence of My on Mps is crucial to 
achieve an accurate determination of the mass of the strange quark. 

In table 2, we present a rather extended set of results for the bare and renormalized 
quark masses, obtained from the vector and axial- vector Ward identities, with the bs com- 
puted in boosted perturbation theory. The renormalized masses, in the RI scheme at /i = 2 
GeV, have been obtained from the bare ones using the renormalization constants of tab. 1. 
We see from this table that, by using quadratic fits, denoted as QQ in the table, and non- 
perturbative Zs, we obtain very consistent results for the two calibrations of the lattice 
spacing, the masses from the vector or axial-vector Wis and the quark mass extracted 
from rriK or m^. In particular, we stress the excellent agreement between the values of 
the quark masses obtained from the vector and axial-vector WI. This agreement is only 
possible if one uses Zs and Zp computed non-perturbatively. This was first noticed in 
ref. 0. A sizeable difference remains, instead, if one uses boosted perturbation theory, see 
also p. For example, with the perturbative values of Zs and Zp, we obtain, at = 2 
GeV, mf^ = 118(12) MeV (from m,) and mf^ = 101(11) MeV (from p,) instead of 135(14) 
MeV and 138(16) MeV, respectively. 

From the results of the table, we extract our best estimate of the light quark masses in 
the RI scheme and at /i = 2 GeV (in the following all the results for the running masses 
refer to = 2 GeV) 

mf^ = 5.6(5) MeV mf^ = 138(15) MeV . (3.15) 
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f3 = 6.2 


Method 




Input 


rfiia 




pia 






(GeV) 






(MeV) 




(MeV) 


p-n L4Q 


2.59(12) 




0.00130(11) 


6.1(3) 


0.00127(11) 


6.0(3) 


K*-K L4Q 


2.69(12) 




0.00120(10) 


5.9(3) 


0.00117(10) 


5.8(3) 


p-TT QQ 


2.84(32) 




0.00108(23) 


5.6(6) 


0.00105(24) 


5.5(6) 


K*-K QQ 


2.83(25) 




0.00109(18) 


5.6(5) 


0.00106(18) 


5.5(5) 


Method 




Input 


rhsa 




Psfl 


nig 




(GeV) 






(MeV) 




(MeV) 


p-TT L4Q 


2.59(12) 


rriK 


0.0312(26) 


147(6) 


0.0313(28) 


149(7) 








0.0363(35) 


171(9) 


0.0379(40) 


181(12) 


K*-K L4Q 


2.69(12) 


rriK 


0.0290(23) 


142(6) 


0.0291(24) 


144(6) 








0.0319(24) 


156(5) 


0.0331(27) 


164(7) 


p-TT QQ 


2.84(32) 


rriK 


0.0261(55) 


135(14) 


0.0261(58) 


137(16) 








0.0260(84) 


134(29) 


0.0265(94) 


139(35) 


QQ 


2.83(25) 


rriK 


0.0263(42) 


135(11) 


0.0263(45) 


137(12) 








0.0264(58) 


136(19) 


0.0270(66) 


141(23) 



Table 2: Improved bare masses of the light (fhia = {fhua + rhda) /2 and pia = {pua + pda)/2) 
and strange (ifisa and pga) quarks in lattice units. The improvement coefficients bm, bA and 
bp necessary to obtain fhi, pi, etc., have been taken from boosted perturbation theory. The 
renormalized quark masses in the RI scheme at a renormalization scale p = 2 GeV, obtained 
using the calibration of the lattice spacing given in this table, and the non-perturbative 
values of the Zs from tab. 1, are also given. We use the values of the inverse lattice 
spacing obtained from mp-m.,^ (denoted as p-Tc) and from mK*-TnK (denoted as K*-K). 
L4Q denotes the same fitting procedure as in ref. J^: a linear fit for My, over the four 
quark masses, combined with a quadratic fit for Mpg. These results can be directly compared 
to those of refs. J^. QQ indicates that quadratic fits were used for both My and Mpg. 
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Note that, to obtain the resuhs in eq. ( p.l5| ), we never used perturbation theory but for 
the values of 6^, &a and bp. However, as discussed above, this is expected to be a source 
of neghgible uncertainty in the determination of hght quark masses. 



4 MS and RGI masses 

Perturbation theory only enters if we want to convert the results of the RI scheme into the 
MS scheme, which has been adopted as the standard one in the literature (both in lattice 
and QCD sum rule calculations). Although this is not necessary (and we believe that m^^^ 
is a more convenient definition, see below), for comparison with other determinations, we 
also give the quark masses in MS. These are found using the relation 

m^(/i) = i?^m^-^(/x) , (4.16) 
where, following the notation of ref. 



R 



with 



(An) V ™yo ' (47r 
« o(^c-l) 



AN, 



c 



(^^Or = ^^96 Np (-75 - 2645iV,^ + 288C3 + 576iV,^C3 + 356iV,n^) (4.18) 

In the above equation, = 1.20206 ... is the Riemann zeta function. Since all previous 

results have been obtained at the NLO, i.e. by ignoring the corrections due to (^Z^^^^^ \ we 

give the results both at the NLO and at the N^LO order. In the numerical evaluation of Rm 
we have used Uf = 4. The reason for this choice, in spite of the quenched approximation 
adopted in our simulation, is the following. The mass m^^{fi) is to be interpreted as 
the mass in the continuum which includes a systematic (unknown) error coming from the 
quenched approximation. Thus, m^\fi) is the estimate of the physical value of the quark 
mass at the scale fi = 2 GeV, at which nf = 4. By fixing in both cases as{Mz) = 0.118, 
which corresponds to af^^{fi = 2 GeV) = 0.296 and af" ^°(/i = 2 GeV) = 0.300, we 
obtain R^^^ = 0.874 and R^^^^ = 0.804 which give 

^NLO MS ^ 4_g(4) j^gy ^ ^NLO MS ^ 121(13) McV 

mf^'^^ = 4.5(4) MeV , mf ^'^ ^ = 111(12) MeV . (4.19) 

Note that by using BPT, we would have been obtained, for example, mf^^ ^'^^ = 4.3 MeV 
and m^^^ '^^^ = 98(11) MeV (from the vector WI, which is the most favorable case for 
BPT). 
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We now compare our results to those of refs. and p. For /i = 2 GeV, we rescale 
their results with the ratios of the perturbative Zs to the non-perturbative ones, at the 
corresponding values of /3. In this way, using the results of ref. 0, we get m^^^ = 126(5) 
MeV from the vector WI and m^^^ = 138(3) MeV from the axial-vector WI and. 



from ref. 



m 



NLO MS 



127(15) MeV from the vector WI. These numbers are in good 



agreement with the NLO result given in eq. ([4.19|) . 

We also analyzed the data of our test run aX (3 = 6.0. In this case we obtain rather larger 
values: m^^'^ ^''^ = 145(7) MeV from the vector WI, using quadratic fits [m^^^ ^'^^ = 
145(14) MeV from the axial WI). With linear fits on the lightest quarks, corresponding to 
K = 0.1324, 0.1333 and 0.1342, we obtain a result much closer to that obtained at f3 = 6.2, 
i.e. m^^^ = 126(6) MeV. As a further check, we have fitted the results of ref. |^ with 
our programs (with quadratic fits), obtaining m^^^ = 131 MeV []. The problem is that 
our accuracy at = 6.0 is rather poor, and we attribute the discrepancy between the 
results with quadratic fits and all the others to the fact that we do not control, in this case, 
the quadratic corrections. These are affected by some differences with the more precise 
results of ref 0, which we have found for the values of the meson masses for the lightest 
quarks. 



A originally discussed in ref. |^ for lattice calculations, and more recently stressed by 
the Alpha collaboration, a more convenient definition of the quark mass is m^^^ which, 
analogously to the RGI 5-parameter, Bk, is also renormalization scheme and scale in- 
variant. Moreover, m^^^ can be directly related to the mass parameters appearing in the 
fundamental Lagrangian of grand unified theories, at the unification scale. The relation 



between m^^^ and m^\n) is given by p9 



RGI 



RRGim^^ (/i) = as (/i) 



-12/25 



as in) /I 7606 \ 
Att V 1875 / 



(4.20) 



al ifi) /3819632767 952 



21093750 



15 



m^^ (/i) 



(4.21) 



With the same values of as given above, we find R'^g? — l-^O and R'^gi'^ — 1-28, corre- 
sponding to 



m 



NLO 

I 



m 



N^LO 

I 



7.8(7) MeV 
7.2(6) MeV 



,NLO 



= 193(21) MeV 
^N^LO ^ 177(19) MeV . 



(4.22) 



5 Conclusions 

Using the O(a^) improved Alpha action, we have computed the light-quark masses renor- 
malized non-perturbatively. The NPM used in this study removes the discrepancies, ob- 
served in previous calculations, between the results obtained using the vector and the 

^ It is not possible to extract an error in this case. 
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axial-vector Ward identities. It also gives values of the quark masses larger than those 
obtained by computing the renormalization constants with (boosted) perturbation theory. 
We use perturbation theory only to translate our results either to the MS scheme, which 
has been adopted in most of the determinations of the quark masses, or to our preferred 
definition which is the RGI mass. Using the results of ref. this can be done at N^LO 
accuracy. 

A further reduction of the systematic errors would require the non-perturbative deter- 
mination of the 6s (although we believe that this uncertainty is, for light quarks, by far 
smaller than the others); the evaluation of Zs{fi) and Zp{fi) at larger physical scales, /i, and 
smaller values of /ia, a fit of the renormalized masses in (in order to remove the residual 
discretization errors) and, of course, precise unquenched calculations. We also found that a 
precise determination of the quadratic dependence of the vector meson mass on the quark 
mass is important to understand the differences in rhg determined using different mesons 
(~ 20 MeV). 



Acknowledgements 

V.L. and G.M. acknowledge the M.U.R.S.T. and the INFN for partial support. We are 
grateful to A. Le Yaouanc for useful and interesting discussions. 



References 

[1] C. R. Allton et al, Nucl. Phys. B431 (1994) 667. 

[2] T. Bhattacharya and R. Gupta, Phys. Rev. D55 (1997) 7203. 

[3] S. Aoki et a/., JLQCD Collaboration, Nucl. Phys. (Proc. Suppl.) 63 (1998) 275; 
CP-PACS Collaboration, Nucl. Phys. (Proc. Suppl.) 60A (1998) 14 and Nucl. Phys. 
(Proc. Suppl.) 63 (1998) 161. 

[4] N. Bicker et al, SESAM Collaboration, Phys. Lett. B407 (1997) 290. 

[5] B.J. Cough et al, Phys. Rev. Lett. 79 (1997) 1622. 

[6] M. Cockeler, R. Horsley, H. Perlt, P. Rakow, C. Schierholz, A. Schiller and P. Stephen- 
son, Phys. Lett. B391 (1997) 388; Phys. Rev. D57 (1998) 5562. 

[7] C. R. AUton, V. Gimenez, L. Ciusti and F. Rapuano, Nucl. Phys. B489 (1997) 427; 
|hep-lat/9801028| . 

[8] A. Cucchieri, M. Masetti, T. Mendes and R. Petronzio, Phys. Lett. B422 (1998) 212. 
[9] K. Jansen et al., Phys. Lett. B372 (1996) 275. 



12 



[10] For a review, see M. Liischer, Summer School in Theoretical Physics, Session 68, Les 
Houches 97, [hep-lat/9802029 . 

K. Symanzik, Nucl. Phys. B226 (1983) 187. 

M. Liischer, S. Sint, R. Sommer and P. Weisz, Nucl. Phys. B478 (1996) 365. 

M. Liischer, S. Sint, R. Sommer, P. Weisz and U. Wolff, Nucl. Phys. B491 (1997) 323. 

M. Liischer, S. Sint, R. Sommer and H. Wittig, Nucl. Phys. B491 (1997) 344. 

M. Guagnelli and R. Sommer, Nucl. Phys. (Proc. Suppl.) 63 (1998) 886. 

G. Martinelli, C. Pittori, C.T. Sachrajda, M. Testa and A. Vladikas, Nucl. Phys. B445 
(1995) 81. 

G. Parisi, in "High Energy Physics 1980", Proceedings of the XXthe International 
Conference, Madison, Wisconsin, eds. L. Durand and L.G. Pondrom (American Insti- 
tute of Physics, New York, 1981). 

G.P. Lepage and P. Mackenzie, Phys. Rev. D48 (1993) 2250. 

D. Becirevic et ai, in preparation. 
M. Bochicchio et ai, Nucl. Phys. B262 (1985) 331. 
L. Maiani and G. Martinelh, Phys. Lett. B178 (1986) 265. 

B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259 (1985) 572. 

G. Martinelh et ai, EDIMBURGH 96/28, unpublished; C. Dawson et ai, Nucl. Phys. 
(Proc. Suppl.) 63 (1998) 877. 

G. Martinelli, G.C. Rossi, C.T. Sachrajda, S. Sharpe, M. Talevi and M. Testa, 
Phys. Lett. B411 (1997) 141. 

G.M. Divitiis and R. Petronzio, Phys. Lett. B419 (1998) 311. 

C. Dawson et ai, in preparation. 

T. Bhattacharya, R. Gupta and W. Lee, presented by W. Lee at the Lattice '98 
conference. 

S. Sint and P. Weisz, Nucl. Phys. B502 (1997) 251. 

E. Franco and V. Lubicz, ROME-I-1 198-98, |liep-ph/980349l| , to appear in Nucl. Phys. 

A. Gonzalez Arroyo, G. Martinelh and F.J Yndurain, Phys. Lett. 117B (1982) 437. 

M. Liischer et al, the ALPHA collaboration, Nucl. Phys. B389 (1993) 247; ibid. B413 
(1994) 481. 



13 



